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ABSTRACT 
A  Markov  type  model  for  studying  rank  distribution  in 
an  hierarchical  organization  is  examined.   Various  sets  of 
rank  distributions  are  defined  and  their  properties  discussed, 
A  computer  aided  test  is  developed  for  testing  a  given  rank 
distribution  to  determine  if  it  is  an  element  of  a  particular 
set.   For  a  distribution  not  in  the  Steady  State  set,  a  test 
is  developed  to  determine  if  the  distribution  can  be  returned 
to  in  m  transitions  or  steps. 
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I.   INTRODUCTION 

Several  authors  have  studied  models  for  analyzing  rank 
distributions  in  an  hierarchical  organization.   Bartholomew 
[1967]  uses  Markov  models  to  a  great  extent  in  his  analysis 
of  various  social  systems.   Thonstad  [1968]  developed  a 
Markov  type  model  for  use  with  the  Norwegian  educational 
system.   More  recently,  Branchflower  [1970]  utilized  a  Markov 
type  model  to  examine  the  promotion  and  appointment  policies 
of  a  university  and  the  resulting  rank  distributions  of 
faculty.   Branchflower  assumed  the  promotion  policies  were 
fixed  and  used  his  model  to  predict  the  faculty  distribution 
in  the  next  time  period  using  the  current  faculty  distributions 
and  actual  appointments  to  each  rank.   In  his  analysis, 
Branchflower  defines  three  sets  of  rank  distributions:  (1)  an 
Attainable  set,   (2)  the  Maintainable  set,  and   (3)  the 
Containment  set.   An  Attainable  set  is  a  set  of  rank  distri- 
butions which  can  be  reached  from  some  previously  attainable 
set.   The  Maintainable  set  is  better  known  as  the  Steady 
State  set  of  rank  distributions.   The  Containment  set  is  the 
set  of  rank  distributions  which  can  be  returned  to  infinitely 
often  by  means  of  a  sequence  of  appointment  policies. 

The  analysis  in  Branchflower  [1970]  leads  us  to  pose 
several  questions  about  the  rank  distributions  given  some 
promotion  policy; 

1.    Can  the  present  rank  distribution  be  maintained  for 
all  time  by  means  of  appointment  policies  alone? 


2.  Can  another  given  rank  distribution  be  attained  from 
the  present  distribution  in  a  finite  time  period,  or  can  the 
desired  distribution  never  be  reached? 

3.  What  are  the  rank  distributions  which  can  be  returned 
to  at  some  time? 

In  this  paper,  our  purpose  is  first,  to  define  explicitly 
the  sets  of  rank  distributions  of  interest  associated  with 
the  model,  and  to  demonstrate  various  properties  of  these 
sets.   Second,  we  shall  use  these  definitions  and  properties 
to  develop  test  procedures  to  determine  if  a  given  distribu- 
tion is  an  element  of  a  particular  set.   This  analysis  will 
provide  answers  to  questions  1  and  2,  and  will  permit  us  to 
give  crude  bounds  as  an  approximation  to  answer  question  3. 
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II.   THE  BASIC  MODEL  AND  NOTATION 
Consider  a  discrete  time  Markov  model  with  n+1  states 
where  state  j  is  the  rank  of  a  member  of  an  hierarchical 
organization.   The  set  of  states  can  be  partitioned  into  two 
subsets:   the  first  n  states  are  transient  states  and  are 
called  active  ranks;  the  second  subset  is  the  n+lst  state 
which  is  an  absorbing  state.   Let  p-  •  be  the  fraction  in 
state  j  at  time  t+1  that  were  in  state  i  at  time  t.   Then 


n+1 

Z   p..  =  1,  1=1,2, ... ,n+l  . 

j=l  13 


We  represent  these  fractions  in  a  transition  matrix 


Q  = 


P    w 
0    1 


where  P  is  an  n  x  n  matrix,  tp^-;]/  i/ j=l/2, . . . ,n,  and  w  is 

an  n  dimensional  column  vector  called  the  loss  vector  with 

w-  =  p •    , .   P  is  called  the  Promotion  matrix  and  describes 

the  transition  process  of  the  active  ranks  in  the  organization 

n 
It  has  the  property  that   E   p-j  <_  1,  i=l,2,...,n  with  strict 

j=l   J 
inequality  holding  for  at  least  one  row. 

We  assume  that  the  organization  is  in  an  equilibrium 
condition,  and  thus  keep  the  total  number  in  the  organization 
constant.   We  normalize  by  assuming  this  number  is  one. 

Let  x(t)  be  an  n  dimensional  row  vector,  x(t)  =  (x, (t), 
x2 (t) , . . . ,xR (t) )  where  x.(t)  is  the  fraction  in  the  ith  rank 
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n 

at  time  t.   Note  that  E   x±(t)  =  1,  x±  (t)  >_  0.   Let  f  (t)  be 

i=l 

an  n  dimensional  row  vector,  f  (t)  =  (f  -.  (t)  ,  .  .  . ,  f  (t) )  where 

f^(t)  is  the  fractional  number  of  new  appointments  to  rank 

i  at  time  t,  f  ^  (t)  _>  0 . 

The  rank  distribution  at  time  t+1  is  dependent  on  the 

rank  distribution  at  time  t,  the  Promotion  matrix,  and  the 

appointment  policy  at  time  t+1.   In  vector  notation, 

x(t+l)  =  x(t)-P  +  f(t+l).  (1) 

The  requirement  that  the  total  number  in  the  organization  be 

constant  requires  that  the  total  fraction  appointed  to  the 

active  ranks  be  equal  to  the  total  fraction  lost  to  the 

n 
absorbing  state:   £   f • (t+1)  =  x(t)*w.  (2) 

i=l 

We  use  the  following  notation  throughout  the  remainder 

of  the  paper:   e  is  an  n  dimensional  vector  with  all  components 

equal  to  1;  e-  is  an  n  dimensional  unit  vector  with  the  jth 

component  equal  to  1 .   It  is  apparent  by  the  context  whether 

they  are  row  or  column  vectors . 
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III.   SET  DEFINITIONS 
The  set  of  all  rank  distributions,  F,  is  defined  as 

F  =  {x;  x*e  =  1,  x  >_  0} 

Note  that  the  time  notation  is  omitted  when  it  does  not  affect 
the  definition. 

An  Attainable  set  is  a  set  of  rank  distributions  which 
can  be  reached  in  one  step  from  a  previous  set  of  rank  distri- 
butions.  Let  m  be  the  number  of  steps  (time  periods)  from 
the  set  of  all  rank  distributions,  F,  then 

A^  =    {x;  X  >  Z'P,  x  e  F,  z  e  Am-l^  '  m  =  i/2****   and 
AQ  =  F. 

Consider  the  sequence  of  Attainable  sets,  {A } .   If  the 
sequence  converges  to  a  set  C  we  call  the  limiting  set  the 
Containment  set.   Mathematically,  we  have 

C  =   lim  Am 

m  ->  °° 

=   lim  (x;  x  >  z'P,  x  e  F,  z  £  ^a-i^ 

m  ■>  °° 

=  {x;  x  _>  z*P,  x  e  F,  z  e   C},   if  the  limit  exists. 

This  definition  is  somewhat  unusual  in  that  the  Containment 
set  is  defined  in  terms  of  itself. 

The  set  of  all  distributions  which  can  be  maintained  in 
successive  time  periods  is  called  the  Steady  State  set,  M. 
Let  f(t)  =  f  for  every  t,  then  since  x(t+l)  =  x(t),   from 
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(1)  we  have  x(t)  =  x(t)«P  +  f 
x(t)  •  (I-P)  =  f  >  0 

so  we  define 

M  =  {x;  x»  (I-P)  >_  0,  X  e  P} 

If  a  rank  distribution  can  be  reached  from  the  Steady 
State  set,  M,  then  we  call  a  set  of  such  distributions  a 
Reachable  set.   All  distributions  which  are  reachable  from  M 
in  m  steps  are  in  the  set  Bm/  where 

Bm  =    {x;  x  j>  z»P,  x  e  F,  z  e  Bm-l^  »  m  =  1/2,...   and 
BQ  =  M. 

Now  consider  the  sequence  of  Reachable  sets.  {B  }.   If 

m 

this  sequence  converges  to  a  set  B,  we  call  the  limiting  set 
the  Reached  set.   Mathematically,  we  have 


B  =   lim  Bm 

m  ■*■  °° 


=   lim  {x;  x  >^  z*P,  x  e  F,  z  e  Bm_-^} 

m  ■*  °° 

=  {x;  x  >  z*P,  x  £  F,  z  e  B}  ,    if  the  limit  exists. 


Thus  by  their  definitions,  we  see  that  the  Containment  set  and 
the  Reached  set  must  be  the  same  set  when  they  both  exist.  In 
section  IV  we  prove  they  always  exist. 
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IV.   SET  PROPERTIES 

In  this  section,  we  list  several  properties  of  the  sets 
defined  in  section  III  without  proofs.   We  then  prove  several 
lemmas  and  two  Theorems  which  establish  the  existence  of  the 
Containment  set  and  Reached  set. 

PROPERTY  1:   The  set  of  all  rank  distributions,  F,  is 
convex  and  closed. 

PROPERTY  2:   Any  Attainable  set,  A^,  is  convex  and  closed. 

PROPERTY  3:   The  Steady  State  set,  M,  is  convex  and 
closed. 

PROPERTY  4:   Any  Reachable  set,  B  ,  is  convex  and  closed. 

J  m 

LEMMA  1:   An  Attainable  set,  A^,  is  a  subset  of  every 
previously  attainable  set. 

Proof:   F  =  AQ  =  {x;  x»e  =  1,  x  _>  0} 

A,  =  {x;  x  >_  z«P,  x  e    F,  z  e  Aq} 
Pick  x  e  A-,  ^  x  £  Aq   therefore  A-,  £  AQ  =  F 
Assume  Am_1  =  A,^  =  ...  =  Ax  c  AQ. 
Define  A  (z)  =  {x;  x  >_  z*P,  x  e   F}  ,  z  e  Am_i   so 

Am  =   U    VZ)  "  Pick  Z*  £  V-l'  then  since  Ani-l  £  V-2' 
z  e  Am-1 

Am(z*)    =    {x;    x  _>    z*«P,    x    e    F,    z*    e   A^_2^    »    and 

Vz*>  =  Vr 


s°  %  =        U        Vz>    =  Ara-1 


z    £    A 

m-1 


LEMMA  2:   The  Steady  State  set,  M,  is  a  subset  of  any 

Attainable  set.  A  . 

'   m 
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Proof:   M  =  {x;  x  _>  x»P,  x  e  F}.   M  £  A  =  F  by  definition. 

Assume  M  -  A   n  .  Pick  z*  e  M  =*  z*  e  Am   -,    and  z*  >  z**P. 
m-1  "i-J.        — 

Therefore,  z*  e  ^  by  definition  of  A^,  and  M  -  Am. 

LEMMA  3:   Any  Reachable  set,  Bm,  contains  every  previously 
reachable  set. 

Proof:   M  =  BQ  =  {x;  x  _>  z»P,  x  e  F>  .   Trivially,  BQ  Q    B-^. 
Assume  that  Bm_2  £  Bm-1  wnere  Bm-1  =  {*/*  x  >_  z*P, 

x  e  F,  z  e  Bm-2-^  ■ 
Define  the  following  sets:  B   ^  (z)  =  {x;  x  >^  z«P, 
x  e  F},  z  e  Bm-2'  an<^  Bm(z)  -  ^x/  x>^z«P,  x  £  F}, 
z  e  %_!.   Pick  z*  e  Bm_2.   Since  Bm_2  £  Bm-1'  z*  e  Bm-1' 
so  that  Bm_^(z*)  =  Bm(z*) .   Since  z*  is  arbitrary,  and 

since  Bm_2  £  Bm-i,  Bm_x  =    U   Bm_1(z)£     U   Bm(z) 

z  e  B^-l        z  e  Bm_1 

=  B  . 
m 

LEMMA  4:   If  w^  >  0  for  at  least  one  i  and  n  >  1,  then 

the  Steady  State  set,  M,  contains  an  infinite  number  of  points. 

Proof: 

(i)   Non-empty.   Assume  that  M  is  empty.   Then  there 

does  not  exist  an  x  such  that  x*  (I-P)  >^  0.   Because 

n 

£   p.  .  <_   1,  i=l,2,...,n  with  strict  inequality  holding 

j=l   D 

for  at  least  one  i,  (I-P)    exists  with  all  elements 

non-negative.   Let  N  =  (I-P)~  .   If  M  is  empty,  then 

there  does  not  exist  an  x  such  that  x  =  f«N  >_  0.   But 

n 
x-  =  E   fj*n. .  and  since  f-  >  0  and  n  •  •  >  0  for  every 

:   i=1  1  id        i  -      id  - 

i/j#  xj  >^  0  and  therefore  M  is  non-empty. 
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(ii)  Infinite  Interior.   Assume  that  M  contains  only 

one  point,  then  x  (t  +  1)  =  x  (t)  for  every  t  and  any 

appointment  vector.   We  rewrite  equation  (1)  as  x(t)  = 

x(t)«P  +  f(t  +  1)  =  f(t  +  l)'(I-P)"1.   We  also  require 

that  f (t  +  1) »e  =  x(t)«w.   We  can  construct  the 

appointment  vector  in  the  following  manner:   let 

f  (t  +  1)  =  cue,  +  ct^e,.  +...+  a  e  where  £a   =  x(t)*w. 
1J-     2  2        nn       ^  i 

Then   f(t  +   1)    =    (oUfOU/ •  •  • /Ot    )    and 

x(t)    =    (a-,,a0,...,a    )  •  (I-P)"1  (3) 

±      2  n 

Consider  the  appointment  vector  f (t  +  1)  when 
a,  =  x(t)»w  and  a^  =  013  =    ...  =  a   =  0 ,  then  x(t)  equals 
x(t)*w  multiplying  the  first  row  of  (I-P)   .   Since 
any  f  which  solves  (3)  generates  a  point  in  M,  consider 
the  appointment  vector  f(t  +  1)  when  012  =  x(t)*w  and 
a^  =  a..  =  ...  =  aR  =  0 ,  then  x(t)  equals  x(t)*w  multiply- 
ing the  second  row  of  (I-P)--'-.   Since  x(t)  is  one  point, 
and  x(t)*w  is  a  scalar,  the  first  two  rows  of  (I-P)"1 
are  identical.   This  implies  that  r((I-P)_1)  <  n  and 
that  ((I-P)""  )    does  not  exist  which  is  a  contradiction 
Thus  since  M  is  non-empty,  M  contains  more  than  one 
point.   Since  M  is  convex  and  closed,  it  contains  an 
infinite  number  of  points. 
THEOREM  1:   The  Containment  Set,  C,  and  the  Reached  Set, 
B,  exist,  are  identical,  and  non-empty. 
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Proof:   Consider  the  sequence  of  sets  {A_} .   By  Lemmas 

1  and  2,  Hc  ...=  L  =A     ^  .  .  .  £  A,  s=  F .   By  Lemma 
'  ^    m-1  1 

4  M  =|=  <j>  =>  A=f<l)rn  =  l/2/.««  •   Thus  we  have  a  mono- 
tone decreasing  sequence  of  sets  which  is  bounded  below, 
therefore  the  limiting  set  does  exist  and  is  non-empty. 
Now  consider  the  sequence  of  sets  {B  }  where  B0  =  M. 

By  Lemma  3 ,  M  £  Bn  £  B-,  Q     .  „•  ;--.c  B„,  ,  c  B   c c  F . 

J         '       u     l  m-1    m 

Thus  we  have  a  monotone  increasing  sequence  of  sets 

bounded  above  by  F.   Therefore  a  limiting  set  exists. 

As  noted  in  section  III  the  Containment  set,  C,  and 

the  Reached  set,  B,  are  the  same  set. 

THEOREM  2:   The  Steady  State  set,  M,  any  Attainable  set, 

A^,  the  Containment  set,  C,  any  Reachable  set,  Bm,  and  the 

set  of  all  rank  distributions,  F,  are  identical  if  and  only 

if  the  Promotion  matrix,  P,  is  diagonal. 

Proof:   Since  M  s  B   £  C  -  K^    e  F/  we  need  only  show  that 

M  =  F. 

(i)    Assume  that  P  is  diagonal.   We  can  write  the 

distribution  vector  as  x(t)  =  a  en  +  a  e   +...+  a  e  , 

1  1    2  2        n  n 

where  a.  >  0.   If  la.  =  1  then  x(t)  e  F.   Rewrite 
i  -        11 

equation  (1)  x(t  +  1)  =  x (t) »P  +  f (t  +  1) .   We  require 
f(t  +  1)  *e  =  x(t)  •  (I-P)«e.   Since  P  is  diagonal,  x(t)« 
(I-P)  is  a  suitable  choice  for  f (t  +  1),  (f . (t  +1)  >  0) . 
So  we  can  now  write  x(t  +  1)  =  x(t)»P  +  x(t)«(I-P)  =  x(t) 
Now  x(t  +  1)  =  x(t)  and  x  (t)  e  F  =»  x  (t)  e  M,   .".   F  £  M. 
And  since  M  e  F,  M  =  F. 
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(ii)    Assume  M=  F.   x(t  +  1)  =  x(t)  for  all  t. 
Assume  that  P  is  not  diagonal.   x(t)  =  x(t)*P  +  f(t  +  1)  . 
since  x(t)  e  F,  we  can  write  x(t)  =  oue.  +  a9e9  +  ...  + 
anen,  a.  >  0,  Za-  =1.   In  particular,  let  on  -  1, 


a^  =  0,  i  =)=  1.   So  x(t)  =  a-j^e-^  and  (1,0,.  ..,0) 


(1,0,. ..,0) 


Pll   Pl2 


'In 


nl   ^n2 


nn 


+   f(t  +  1) 


-  (P- 


plr.)  +  (f-i  ,f9  r.../  f  )•   This  requires 


'11'  ^12  '*•*'  ^ln'    VJ-1/J"2        n- 
that  the  following  equations  be  satisfied: 


fl   +   Pll   -   1 


2      +      Pl2 


=   0 


n 


Pin 


=   0 


Since  f  .  ^  0 ,  p.  .  ^  0,  we  have  f.  =  p. .  =  0,  j  =f  1.   We 
repeat  this  procedure  for  a-  =1,  ak  =  0,  k  =|=  j  ,  j  = 
2,3,...,n.   Thus  p^    =  0,  i  =}=  j  and  P  is  diagonal. 
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V.   EXTREME  POINTS  OF  THE  VARIOUS  SETS 
The  characterizations  of  the  various  sets  in  the  model 
in  the  above  section  do  not  describe  them  in  a  useful  manner. 
If  we  knew  the  extreme  points  of  all  of  the  sets,  then  we 
would  have  almost  all  of  the  information  about  the  sets, 
since  any  interior  point  of  a  set  is  simply  the  convex  combin- 
ation of  its  extreme  points. 

For  two  sets  in  particular,  the  extreme  points  are  easily 
obtained.   The  set  of  all  rank  distributions,  F,  is  the 
fundamental  simplex  in  En  with  extreme  points,  e^,  i  =  1,2, ...n, 
The  subset  of  F  which  is  most  easily  described  is  the  Steady 
State  set,  M.   Each  vector  in  the  set  satisfies  the  inequality 
x  2l  x#P.   This  can  be  written  in  equation  form  as  x*  (I-P)  =  f 
where  f  is  the  vector  of  appointments  to  the  different  ranks. 
We  define  N  =  (I-P) "^  so  we  write  x  =  f*N.   Since  we  require 

x*e  =  1  and  since  f  can  be  written  f  =  Ed£  (x«  (I-P)  »e)  ,  aj_  >_  0, 

i 

Ea^  =  1,  we  write  the  equation  in  an  equivalent  form  as  x  =  f*B 
i 

Nii 
where  f  =  Za-e.  and  B •  .  =  - — +-   .   That  is,  the  extreme  points  of 

J  3 

the  Steady  State  set  are  given  by  the  rows  of  (I-P)    after  the 
rows  are  normalized. 

The  problem  of  determining  the  extreme  points  of  the  other 
sets  in  the  model  is  much  more  difficult.  If  we  are  working  in 
three  dimensions,  it  is  not  too  difficult  to  calculate  all  of 

the  extreme  points  of  one  of  the  sets,  A  or  B  ,  for  a  small  m. 

c  mm 
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Each  set  in  the  model  has  the  property  that  each  point  in  the 
set  generates  a  subsimplex  in  the  fundamental  simplex.   The 
extreme  points  of  the  generated  set  are  found  from  the  convex 
hull  of  the  subsimplices  generated  by  the  extreme  points  of 
the  original  set.   In  three  dimensions  we  use  a  geometrical 
method  to  determine  the  extreme  points  of  the  convex  hull  of 
the  generated  subsimplices.   If  we  know  the  extreme  points  of 
^  for  some  m,  then  we  project  Am  onto  a  plane  generated  by 
two  unit  vectors.   Then  we  calculate  the  extreme  points  of 
the  subsimplex  generated  by  each  extreme  point  of  Am  and  project 
these  onto  the  plane.   Then  we  simply  test  linear  combinations 
of  the  generated  extreme  points  to  determine  that  there  are 
no  points  which  lie  beyond  the  linear  combination  in  each  of 
the  unit  directions.   The  points  which  satisfy  this  criterion 
are  the  extreme  points  of  the  set  Am+i .   A  similar  procedure  is 
used  for  Bm. 

When  P  is  upper  triangular,  some  of  the  extreme  points  of 
AJt^  are  also  extreme  points  of  Am+^  (e.g.  en  is  an  extreme 
point  of  Am  for  every  m.)   These  remain  extreme  points  for  all 
succeeding  sets  generated  from  Am+2.   Similarly,  Bm  has 
extreme  points  common  to  all  ^  which  are  extreme  points  of 
A  as  well,  and  therefore  are  extreme  points  of  the  Containment 
set,  C.   Figure  1   depicts  several  attainable  sets  drawn  on 
barycentric  coordinates  for  a  three  dimensional  rank 
distribution  where 
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FIGURE  1.   ATTAINABLE  SETS  ON  BARYCENTRIC  COORDINATES 
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p  = 


.5  .4  0 
0  .6  .3 
0      0.7 


Table  I   contains  the  coordinates  of  the  numbered  points  on 
Figure  1.   The  heavy  line  in  the  figure  represents  the  known 
portion  of  the  boundary  of  the  Containment  set  after  calculating 
A7 .   The  portion  of  the  boundary  between  points  9  and  10 
changed  very  little  between  Ag  and  A7  so  the  boundary  of  Ay 
is  probably  a  good  approximation  of  the  boundary  of  the 
Containment  set  in  this  case. 

TABLE  I:   POINTS  ASSOCIATED  WITH  FIGURE  1. 


Point 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Rank 

1 

0 

.3 

.39 

.393 

.3723 

.35289 

.18615 

.176445 

.0882225 

0 

1/3 

2 

0 

0 

.12 

.228 

.2940 

.32532 

.49206 

.500706 

.5355714 

.5 

1/3 

3 

1 

.7 

.49 

.379 

.3337 

.32179 

.32179 

.322849 

.3762061 

.5 

1/3 

It  is  not  clear  how  this  graphical  projection  method  can 
be  extended  to  higher  dimensions  to  obtain  the  extreme  points. 
Nor  is  it  clear ,  based  on  the  above  example,  whether  or  not 
the  Containment  set  has  a  finite  number  of  extreme  points. 
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VI.   THE  DEVELOPMENT  OF  SEVERAL  TESTS 
Since  it  is  not  feasible  to  analytically  determine  the 
extreme  points  of  all  of  the  sets,  we  must  devise  some  method 
to  answer  our  original  questions.   The  answers  to  those 
questions  depend  on  knowing  the  extreme  points  of  all  of  the 
sets.   To  provide  partial  answers  to  the  questions,  we  develop 
test  procedures  to  determine  if  a  given  vector  is  an  element 
of  a  particular  set,  and  to  determine  if  a  given  vector  can 
be  returned  to  in  m  steps  or  reached  from  another  specified 
vector  in  m  steps.   We  know  that  a  given  vector  cannot  be 
returned  to  at  any  time  if  there  exists  an  i  such  that  the 
test  vector  is  not  in  A^. 

Each  transition  from  one  set  to  another  requires  that  a 
certain  set  of  inequalities  be  satisfied.   When  we  consider 
the  system  of  inequalities  that  account  for  all  of  the  transi- 
tions, we  need  only  test  for  feasibility.   If  no  feasible 
solution  exists,  then  the  test  vector  fails  the  test.   Several 
specific  tests  are  developed  in  the  next  sections. 


A.    A  TEST  FOR  Bm 

m 

In  this  section,  we  develop  a  test  to  determine  if  a  given 
rank  distribution  is  an  element  of  the  set  Bm.   Recall  that 

Bm  =  {zm;  zm  1   zm-l#p'  zm  e  F'  zm-l  e  Bm-1}  for  m  =  1>2"- 
and  that  B   is  the  set  of  rank  distributions  which  is  reachable 

in  m  steps  from  the  Steady  State  set.   Row  vectors  have  been 

used  in  all  developments.   To  simplify  notation,  the  test 

inequalities  will  be  specified  in  tableau  form  and  vectors 
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will  be  column  vectors.   Transposition  is  indicated  for  the 
prombtion  matrix  where  it  appears . 

We  now  write  the  sequence  of  inequalities  which  must  be 
satisfied  if  the  vector  X  is  in  the  set  Bm  as 


z0     zl    z2   *   *   *    Zm-2     Zm-1 
(PT-I)  <  0 


PT       -I  <  0 


PT    -I  <  0 


PT        -I     £  0 
PT      <  XT 


(4) 


e  =1 

e  =1 

e  =1 

e      =1 

The  constraint  matrix  appears  deceptively  smaller  than  it 
is.   Since  each  z.^   is  an  n  vector,  the  dimensions  of  the  matrix 
are  mn  +  m  +  n  rows  and  mn  columns .   Even  though  the  constraint 
matrix  is  very  sparse,  we  were  unable  to  apply  any  algorithm 
which  would  process  the  problem  more  efficiently  than  the 
Simplex  algorithm. 
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It  is  clear  that  a  feasible  solution  to  (4)  is  simply  a 
Phas'e  I  Simplex  solution.   The  method  which  we  use  to  obtain 
this  solution  is  to  add  an  artificial  variable  to  each  of  the 
feasibility  constraints,  the  last  m  rows.   We  then  minimize 
the  sum  of  the  artificial  variables.   A  feasible  solution  to 
the  new  system,  (4)  plus  the  artificial  variables,  always 
exists.   The  result  of  the  test  is  determined  by  the  values  of 
the  artificial  variables.   If  all  of  the  artificial  variables 
are  equal  to  zero,  then  (4)  has  a  feasible  solution  and  we 
conclude  that  the  test  vector,  X,  is  an  element  of  Bm.   If 
any  of  the  artificial  variables  are  not  zero,  then  a  feasible 
solution  to  (4)  does  not  exist  and  we  conclude  that  X  is  not 

in  Bm* 

B.    A  TEST  FOR  P^ 

In  this  section,  we  develop  a  test  to  determine  if  a 
given  rank  distribution  is  an  element  of  Am,  i.e.  if  it  can 
be  reached  in  m  steps  from  some  point  in  F.   Recall  that 

^n  =  {zm'-  zm  1  zm-l'P'  zm  e  F'  zm-l  e  V-l}  for  m  =  1'2"" 
We  write  the  set  of  inequalities  as 
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z0      zl     z2   *   '   *    zm-2      zm-l 


-I  <  0 


PT     -I  <  0 


PT        -I      <  0 

pT      <  XT   (5) 

e  =1 

e  =1 

e  =1 

e      =1 

The  constraint  matrix  has  ran  +  m  rows  and  mn  columns.   We 
now  test  (5)  for  a  feasible  solution  in  the  same  manner  as 
the  test  for  Bm.   The  results  of  the  test  are  interpreted  in 
the  same  manner  as  those  for  Bm.   Note  that  the  constraint 
matrix  in  (5)  is  the  same  as  that  in  (4)  without  the  first  n 
rows . 

C.    A  TEST  FOR  RETURN  IN  m  STEPS 

In  this  section,  we  develop  a  test  to  determine  whether 
or  not  a  given  rank  distribution,  X,  can  be  returned  to  in  m 
steps.   If  X  is  an  element  of  the  Steady  State  set,  we  know 
that  it  can  be  returned  to  in  one  step,  i.e.  it  can  be 
maintained  forever.   So  we  are  interested  in  a  test  for  m 
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greater  than  one.   Since  there  are  m  transitions,  there  must 
be  m-1  intermediate  vectors  because  the  initial  and  terminal 
vectors  are  both  X.   The  first  intermediate  vector,  z-.,  must 
satisfy  z-,  _>  X»P.   The  second  intermediate  vector,  z2,  must 
satisfy  z2  ^L  z-,«P.   This  process  continues  until  the  last 
intermediate  vector,  zm_i/  must  satisfy  zm, *P  £  X.   We  write 
the  sequence  of  inequalities  in  tableau  form  as 


1       2  m-2      m-1 


>  PT«XT 


PT      -I  <  0 


>T 


<  0 


PT       -I      <  0 


(6) 


PT      <  xT 


=  1 


=  1 


e  =1 

e      =1 

Again  we  notice  that  there  is  a  great  similarity  between 
the  form  of  this  tableau  and  those  for  the  tests  for  K^   and 
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B  .   The  first  n  inequalities  in  (6)  are  different  from  those 
in  (4) ,  as  well  as  there  being  one  less  intermediate  vector. 
The  constraint  matrix  in  (6)  has  mn  +  m  +  n  rows  and  mn  -  n 
columns . 

As  with  the  tests  for  A^  and  B  ,  any  method  of  finding 
a  feasible  solution  to  (6)  is  acceptable.   One  method  would 
be  to  add  artificial  variables  not  only  to  the  last  m  rows 
of  (6),  but  also  to  the  first  n  rows,  then  minimizing  the 
sum  of  the  artificial  variables.   If  all  of  the  artificial 
variables  are  zero,  then  we  conclude  that  the  rank  distri- 
bution, X,  can  be  returned  to  in  m  steps.   The  method  used 
in  section  VII  involved  adding  artificial  variables  only  to 
the  last  m  rows  of  (6) ,  then  minimizing  the  sum  of  the  artifi- 
cial variables.   This  system  may  not  always  have  a  feasible 
solution.   The  test  requires  that  if  the  system  has  a  feasible 
solution,  then  if  all  of  the  artificial  variables  are  zero, 
we  conclude  that  X  can  be  returned  to  in  m  steps. 

D.    A  TEST  FOR  CHANGING  DISTRIBUTIONS  IN  m  STEPS 

We  now  develop  a  test  to  determine,  given  some  distribu- 
tion vector,  X-,  ,  whether  or  not  we  can  reach  another  specified 
distribution  vector,  Y^2>    in  m  steps.   This  is  nearly  identical 
to  the  problem  discussed  in  the  above  section,  i.e.  we  have  a 
fixed  initial  vector  and  a  fixed  terminal  vector.   The  con- 
straint matrix  is  the  same  as  (6)  but  now  the  right  hand  side  is 
[PT«X^,  0,...,  0,  x"?,  1,...,  1]  .   The  test  procedure  is  then 
the  same  as  the  test  for  return  in  m  steps. 
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E.    A  TEST  FOR  EXTREME  POINTS  OF  THE  CONTAINMENT  SET 

The  similarity  between  the  tests  for  A^  and  B   noted  in 
sections  A  and  B  above  permits  us  to  develop  a  procedure  with 
which  we  can  determine  extreme  points  of  the  Containment  set, 

C. 

T 
Let  G^  =  [  (I-P  ) ,  0,...,  0]  be  an  n  by  mn  matrix  and  let 


G2  = 


>T 


-I 


.T 


e 
•e 


e 
-e 


-I 


-I 


e 
-e 


e 

-e 


m 

be  an  mn  +  2m  by  mn  matrix.   Now  let  X_  =  [0, . . . ,0,X  ,1, . . . ,1] 

T       T 
be  an  mn  +  2m  by  1  column  vector  and  Z  =  [z/.../z±   ]  be  an 

0      m-1 
mn  by  1  column  vector.   Then  the  test  for  A^  reduces  to  finding 


30 


a  feasible  solution  to 


2-2 


(7) 


If  there  is  a  Z*  which  satisfies  (7)  then  the  test  for  B 
would  seem  to  result  in  testing  to  see  if  G-,*Z*  >^  0.   However, 
there  is  no  guarantee  that  Z*  is  unique  and  there  could 
exist  a  Z'  which  solves  (7)  and  also  solves  G-^'Z'  >_  0  for  the 
same  test  vector  X.   So  this  would  be  a  sufficient  test  for 
inclusion  in  Bm/  but  not  necessary.   The  only  necessary  and 
sufficient  test  is  that  described  in  section  III. A. 

Now  let  H>L  =  [ (I-PT) ,0, . . . , 0]  be  an  n  by  mn  +  n  matri: 
and  let 


.x 


H2  = 


PT    _j 


e 
■e 


pt    -i 


e 
-e 


-I 

T.T 


-I 


e 
•e 


e 

-e 
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be  an  mn  +  n  +  2m  by  mn  +  n  matrix.   Let  Z'  =  [zn/....z  ]  be 

u      m 

an  mn  +  n  by  1  column  vector  and  let  XA  =  [0,...,0,1,...,1] 
be  an  mn  +  n  +  2m  by  1  column  vector.   This  amounts  to  letting 

the  test  vector.  X,  become  an  unknown,  zr.      Then  all  of  the 

'  m 

basic  feasible  solutions  to  H2»Z'  _<  XL    are  extreme  points  of 
A  .   If  a  point  is  an  extreme  point  of  P^   and  B   simultaneously, 
then  it  is  an  extreme  point  of  the  Containment  set,  C.   Thus 
we  need  only  test  each  basic  feasible  solution  to  see  if 
Hl*z'  ii  0.   If  it  is,  then  Z1  is  an  extreme  point  of  C. 

Due  to  the  computational  difficulties  associated  with 
this  test,  it  is  presented  for  theoretical  interest  only. 
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VII.   A  CASE  STUDY 

The  tests  developed  in  section  VI  above  were  applied  to 
faculty  distributions  in  the  College  of  Engineering  at  the 
University  of  California  at  Berkeley  based  on  data  from  1960 
to  1968.   The  bulk  of  the  data  is  in  Branchf lower  [1970]  and 
only  that  which  is  relevant  is  reproduced  here. 

We  used  the  Mathematical  Programming  System  on  the  IBM 
360/6  7  computer  at  the  W.  R.  CHURCH  Computer  Center,  Naval 
Postgraduate  School  to  conduct  most  of  the  tests.   This  is 
really  a  trivial  exercise  once  the  data  cards  are  prepared 
for  the  control  programs.   One  FORTRAN  IV  program  was  written 
to  generate  the  data  cards  for  the  test  for  Bm.   The  same 
cards  were  used  to  test  for  Am  after  removing  the  cards  per- 
taining to  the  first  n  rows  of  the  constraint  matrix.   Another 
FORTRAN  IV  program  was  written  which  was  used  to  generate 
data  cards  for  testing  return  in  m  steps.   Again,  hand  modifi- 
cation of  the  right  hand  side  permitted  the  same  cards  to  be 
used  to  test  for  changing  distributions  in  m  steps. 

One  other  FORTRAN  IV  program  was  written  which  was  used 
to  determine  the  extreme  points  of  the  steady  state  set.   It 
was  also  used  to  test  to  determine  if  a  given  vector  was  an 
element  of  the  Steady  State  set. 

All  of  the  tests  were  verified  for  accuracy  by  using  a 
test  case,  namely  the  three  state  promotion  matrix  in  section 
V  and  points  which  were  in  known  sets  from  Figure  1. 
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10      11      12      13 


STATE 

IN*      123456789 

1  .28    .56    .04    .04 

2  .41    .47    .07 

3  .44    .37    .14 

4  .38    .52    .04 

5  .47    .43    .10 

6  .57    .41    .01 

7  .64    .34 

8  .67    .30    .01 

9  .69    .25    .02 

10  .71    .24    .01 

11  .73    .25    .01 

12  .77    .20 

13  .97 
FIGURE  2.   PROMOTION  MATRIX  FROM  1960-1968  DATA 


TABLE  II.   TEST  VECTORS 


^X.VECTOR 

stateX. 

Xi(1965) 

X2(1966) 

X3(1967) 

X4(1968) 

1 

.0161427 

.0154639 

.0 

.0047619 

2 

.0802139 

.0567010 

.0628019 

.0190476 

3 

.0534759 

.0979381 

.0821256 

.0809523 

4 

.0374331 

.0257731 

.0628019 

.0761904 

5 

.0534759 

.0618556 

.0483091 

.0523809 

6 

.0481283 

.0463917 

.0628019 

.0666666 

7 

.1390374 

.0979381 

.0821256 

.0952380 

8 

.1176470 

.1494845 

.1400966 

.1047619 

9 

.1122994 

.1134020 

.1159420 

.1238095 

10 

.1390374 

.1030927 

.1062801 

.0904761 

11 

.0962566 

.1185567 

.1014492 

.1047619 

12 

.0855614 

.0670103 

.0676328 

.1047619 

13 

.0213903 

.0463917 

.0717328 

.0761904 
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For  the  case  study,  we  used  the  thirteen  active  states 
of  the  faculty  distribution  and  the  1960-1968  promotion 
matrix  from  Branchflower  [1970].   This  is  reproduced  as 
Figure  2.   The  actual  distributions  of  faculty  for  1965-196! 
were  normalized  and  used  as  the  four  test  vectors.   These 
are  included  in  Table  II. 

The  test  vectors  were  tested  for  being  elements  of  A_ 

3  m 

and  Bm/  m  =  1,2,. ..,7,  and  were  tested  to  determine  if  they 
could  be  returned  to  in  m  steps  for  m  =  2,3,. ..,8.   The 
resulting  location  of  the  vectors  of  rank  distributions  on 
the  fundamental  simplex  is  indicated  in  Table  III. 

TABLE  III.   TEST  RESULTS 


xl 

x2 

x3 

X4 

Xl  e  A2 
Xl  *  A3 

X2  *  Al 

X3  e  A2 
X3  £  A3 

X4  e  A2 
x4  i  A3 

Recall  that 


if  X.  e  A.  then  after  one  step,  X.  ,  >  X- -P, 


X^+2  is  at  least  an  element  of  A.+^.   The  results  in  Table  III 
are  inconsistent  with  this  theoretical  result  and  require  an 
explanation.   There  are  two  reasons  why  this  could  have  occured, 
The  first  reason  is  that  the  Promotion  matrix  was  an  average 
over  an  eight  year  period  and  did  not  exactly  represent  the 
true  Promotion  matrix  for  each  year.   The  second  reason  is  that 
the  faculty  size  did  not  remain  constant,  but  in  fact  was 
increasing.   An  expanding  organization  can  appoint  personnel 
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to  replace  the  losses  and  then  hire  the  additional  number 
for  expansion  into  any  state.   When  P  is  upper  triangular, 
the  extreme  points  of  successive  A  may  tend  to  move  away 
from  the  lower  ranks.   For  a  three  state  example,  consider 
an  initial  rank  distribution  of  (1,  0,  0)  with 


.5 

.4 

0 

0 

.6 

.3 

0 

0 

.7 

then  after  one  transition,  the  rank  distribution  before  new 
appointments  is  (.5,  .4,  0).   Say  the  number  of  persons  in 
the  organization  is  constant  at  100,  then  the  distribution  of 
persons  in  the  ranks  is  (50,  40,  0)  with  a  loss  of  10  to 
replace.   Appointments  could  result  in  the  distributions 
(60,  40,  0),  (50,  50,  0),  or  (50,  40,  10)  or  any  linear 
combination  of  them.   If  the  organization  expands  by  10  persons, 
then  the  distributions  could  be  (70,  40,  0),  (50,  60,  0),  or 
(50,  40,  20)  or  any  linear  combination  of  them.   If  we 
normalize,  then  these  expanded  distributions  are  outside  A-^ 
for  the  constant  size  case.   The  result  is  that  the  possible 
distributions  for  the  expanding  size  model  can  be  outside  of 
the  normal  set  of  distributions  obtained  for  the  constant  size 
model. 

We  then  aggregated  several  ranks  in  the  model  to  determine 
if  the  aggregated  model  behaved  in  the  same  way.   The  first 
four  ranks  are  Assistant  Professors,  the  next  three  ranks  are 
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Associate  Professors,  and  the  last  six  ranks  are  Full 
Professors.   The  resulting  Promotion  matrix  is 


P   = 


.7649 

.1752 

0 

0 

.8030 

.1841 

0 

0 

.9707 

The  normalized  rank  distributions  for  1965-1968  are 
listed  in  Table  IV. 

TABLE  IV.   TEST  VECTORS:   AGGREGATED  MODEL 


\yECTOR 

rmk\ 

X-^1965) 

X2(1966) 

X3(1967) 

X4(1968) 

Asst. 
Prof. 

Assoc. 
Prof. 

Full 
Prof. 

.187166 
.240642 
.572192 

.195876 
.206186 
.597938 

.207729 
.193237 
.599034 

.180952 
.214286 
.604762 

These  vectors  were  tested  to  determine  their  location  on 
the  fundamental  simplex.   The  results  are  shown  in  Table  V. 

TABLE  V.   TEST  RESULTS:   AGGREGATED  MODEL 


*1 

x2 

x3 

x4 

Xl  £  A9 
Xl  *  A10 

X   e  A 
2     7 

X2  *  A8 

X   e  A 
3     6 

X3  *  A7 

X  e  A 
4     9 

x4  r   A10 
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These  results  using  the  aggregated  model  are  similar  to 
the  results  of  the  thirteen  state  model.   We  noted  above  that 
as  faculty  size  increases,  the  rank  distribution  may  move  to 
a  superset  of  the  set  containing  the  original  distribution. 
In  fact,  this  happens  in  1966  for  an  increase  of  seven  faculty 
members  (3.7%),  and  also  in  1967  for  an  increase  of  thirteen 
faculty  members.   In  1968,  there  was  an  increase  of  three 
faculty  members,  but  the  rank  distribution  moved  to  a  subset 
of  the  set  containing  the  1967  rank  distribution. 
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VIII.   CONCLUSIONS  AND  RECOMMENDATIONS 
We  have  established  several  properties  of  several  sets 
of  rank  distribution  in  an  hierarchical  organization  which 
have  not  been  previously  examined.   In  addition,  based  on 
these  properties,  we  developed  several  tests  to  determine  if 
a  specified  rank  distribution  is  an  element  of  certain  sets. 
The  tests  were  applied  to  data  from  a  case  study  of  the 
College  of  Engineering  faculty  at  the  University  of  California, 
Berkeley. 

Several  questions  remain  unanswered  and  are  possible 
areas  for  further  study.   Is  there  some  neighborhood  of  a 
given  rank  distribution  which  can  be  maintained?   Is  there 
an  efficient  algorithm  which  will  solve  the  test  inequalities 
other  than  the  Simplex  method?   Is  there  a  method  of  reducing 
the  number  of  basic  feasible  solutions  which  must  be  found  to 
determine  the  extreme  points  of  the  Containment  set?   Does 
the  Containment  set  have  a  finite  number  of  extreme  points? 
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#  r   +   *   ^   A  4   *  ^(   )((   ^c  ^   *':   H-   A  *  A   *   *  ^   y,:   #   *   s}:   ^   !k   iV  ^:   t   * 

*  *  *  *  DATA  GENERATOR  -  TEST  FOR  RETURN  IN  M  STEPS  *  *  *  * 
******  ROBERT  Lo  ARMACOST  -  -  -  1  JULY  1970  *  *  *  *  *  * 
****************************** 

****************************** 

THIS  PROGRAM  WILL  PUNCH  OUT  A  SET  OF  DATA  CARDS  FOR  A 
LINEAR  PROGRAM  USING  THE  MPS  PACKAGE  ON  THE  IBM  360/67^ 
THE  BASIC  INPUT  TO  THE  PROGRAM  IS  A  SQUARE  PROMOTION 
MATRIX  OF  SIZE  N,  AND  K  VECTORS  WHICH  ARE  TO  BE  TESTED 
TO  DETERMINE  IF  THE  VECTOR  CAN  BE  RETURNED  TO  IN  M 
STEPSo   THE  PROGRAM  WILL  GENERATE  DATA  DECKS  FOR  A 
SERIES  CF  TESTS  TO  RETURN  IN  M  STEPS-,   MSTART  IS  THE 
STARTING  INDEX  FOR  THE  NUMBER  OF  STEPS  TO  RETURN*  AND 
MSTOP  IS  THE  ENDING  INDEXo   THESE  INDICES  ARE  READ  IN 
AS  DATAo   THE  FIRST  VECTOR  TO  BE  TESTED  IS  CALLED 
"TEST"o   EACH  SUCCEEDING  VECTOR  IS  CALLED  "REVI"  WHERE 
I  IS  THE  NUMBER  OF  THE  REVISION  OF  THE  RIGHT  HAND  SIDEo 
THE  VARIABLE  IR  IS  THE  NUMBER  OF  REVISIONS  IN  THE 
CONTROL  PROGRAMo   THE  ARRAYS  SHOULD  BE  DIMENSIONED  IN 
THE  FOLLOWING  WAY  - 

AX( (MSTOP-1 )*(N+1)  +  N,  (MST0P-1)*(N+1 H 

KEY  -  SAME  AS  AX 

X(N, IR  +  1  ) 

P(N,N) 

PT(N,N) 

PX(NtN) 

UNI(N,N)o 

THE  PROGRAM  WILL  ACCOMODATE  A  PROMOTION  MATRIX  WITH 
N  =  99  AND  WITH  MSTOP  <  100o 

****************************** 

****************************** 

THERE  IS  A  CERTAIN  AMOUNT  OF  REPETITION  BETWEEN  THIS 
PROGRAM  AND  THE  FOLLOWING  PROGRAM  FOR  A  TEST  FOR  A 
REACHABLE  SET  IN  M  STEPSo   TO  MINIMIZE  DOCUMENTATION, 
CERTAIN  MODULES  IN  THIS  PROGRAM  NEED  ONLY  BE  REPLACED 
BY  OTHER  MODULES  TO  OBTAIN  THE  NEW  PROGRAMo   THESE 
MODULES  ARE  MARKED  WITH  A  NUMBER  IN  THE  LEFT  HAND 
COLUMN  WHICH  MEANS  THEY  ARE  TO  BE  REPLACED  BY  THE 
STATEMENTS  WITH  THE  CORRESPONDING  NUMBER  IN  THE  NEXT 
SECTION  TO  OBTAIN  THE  NEW  PR0GRAMo   IN  ADDITION,  THERE 
ARE  SEVERAL  STATEMENTS  WHICH  NEED  ONLY  BE  REMOVEDo 
THESE  ARE  MARKED  WITH  AN  R  IN  THE  LEFT  HAND  COLUMNo 

****************************** 

INTEGEP*2    KEY(111,98) 

DIMENSION    AXUllt  98)  ,     X(13,4),     P<13,13),    PT(13,13), 
1UNU13.13)  ,    PX(13,13) 

READ    INPUT    DATA 

1  READ<5,600)    MSTART,     MSTOP,    N,     IR 

600       F0RMATKI5) 
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601 

R 
R 

IRl    =    IR    +    1 

READ(5,601)     (  ( P( I, J),J=1,N) 

FORMAK8FIO0O) 

READ(5,601)     ( (X( I, J) ,1=1, N) 

MBEG    =    MSTAPT    -    1 

MEND    =    MSTGP    -    1 

DO    999    M    =    MBEG, MEND 

INITIALIZE 

1=1, N) 
J  =  l,  IRl) 


2 
2 
2 
2 
2 
2 
2 
2 
2 

R 
R 
R 
R 
R 
R 
R 
R 
R 

3 
3 
3 
3 
3 
3 
3 


ITEST  =  1 
JLIM  =  M*N  +  M 
ILIM  =  JLIK  +  N 
DO  10  J=1,JLIM 
DO  10  I  =1, ILIM 
AX (I  ,J)  =  Oo 
10    KEY( I, J)  =  C 

CALCULATE  THE  TRANSPOSE  OF  THE  PROMOTION  MATRIX,  PT , 
AND  LOAD  THE  IDENTITY  MATRIX  UNI 

DO  20  1=1, N 
DO  20  J=1,N 
UNK  I,  J)  =  Co 
PT(I  ,J)  =  P( J, I) 
IF(IoEOoJ)  UNI(I,J)  =  lo 
20  CONTINUE 

CALCULATE  THE  RHS  FOR  THE  FIRST  N  ROWS 

DO  1  J  =  1,IR1 

DO  2  I  =  1,N 

SUM  =  Co 

DO  3  K  =  1,N 
3  SUM  =  SUM  +  PT(I,K)  *  X(K,J) 
2  PXU  ,J)  =  SUM 
1  CONTINUE 


INSERT  UNI  IN  THE  FIRST  N  ROWS  AND  COLUMNS  OF  THE 
CONSTRAINT  MATRIX,  AX 

DO  30  1=1, N 
DO  30  J=1,N 
AXU  ,J)  =  UNK  I, J) 
30  IF(AX(I,J»oNEoO)  KEYU, J)  ■  1 
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50 


INSERT  PT  AND  -I  IN 
(M-2)*N  ROWS  OF  AX 

LASTC  =  N 

NXTR  =  N+l 

NXTC  =  N+l 

JLASTC  =  1 

LASTR  =  2*N 

JNXTC  =  2*N 

MOD  =  1 

IF(MODoEQo 

DO  50  I  = 

DO  40  J  = 

AXU, J)  = 

IF(AX( I , J) 

DO  50  J  = 

IF(I  0EO0J)  AXU  ,J) 

IF(AX( I , J)oEQo-lo0) 

NXTC  =  NXTC  +  N 

NXTR  =  NXTR  +  N 

LASTC  =  JNXTC 

JLASTC  =  NXTC  -  N 

LASTR  =  NXTR  +  N  - 

JNXTC  =  LASTC  +  N 

MOD  =  MOD  +  1 

GO  TO  5 


THE  APPROPRIATE  COLUMNS  OF  THE  NEXT 


M)  GO  TO  100 

NXTR, LASTR 

JLASTC, LASTC 

PT( I-M0D*N,J-(M0D-1)*N) 

0NE0O)  KEYU  ,J)  =  1 

NXTC, JNXTC 

=  -  loO 
KEYU  ,J)  =  1 
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INSERT  PT 
AX 


IN  THE  CORRECT  COLUMNS  OF  THE  NEXT  N  ROWS  OF 


100 


60 


65 

70 

700 
701 
702 

703 

350 
753 
200 


351 
741 

352 
742 

353 
743 

354 
744 
201 


705 

320 
720 
202 


00  60  I  =  NXTRtLASTR 

DO  60  J  =  JLASTCtLASTC 

AX(I,J)  =  PT( I-MOD*N, J-(M0D-1 >*N) 

IF( AX( I , J)oNEo0)  KEY(ItJ)  =  1 

INSERT  l'S  IN  THE  CORRECT  COLUMNS  OF  THF  LAST  M  ROWS 
AX  WHICH  CORRESPOND  TO  THE  FEASIBILITY  CONSTRAINT 

MOD  =  MOD  +  1 
IM  =  LASTR  +  1 
IMM  =  LASTR  +  M 
DO  70  I=IM,IMM 
ICOL  =  (  I-IM)*N  +  1 
ICOL  -  1  +  N 
ICOL,JCOL 


OF 


JCOL  = 
DO  65  J 
AX( I  ,J) 
KEY(  I, J 
K  =  M*N 
AX( I ,K) 
KEY( I,K) 
CONTINUE 


)    = 


lo 

1 

I  - 

lo 

1 


LASTR 


PUNCH    OUT    THE    DATA    SET    NAME    AND    ROWS    CONSTRAINTS    CARDS 


WRITE(6, 
FORMAT( • 
WRITE16, 
FORMAT ( • 
WRITF(6, 
FOPMAT(T 
DO  200  J 
IFUoGTo 
WRITE(6, 
FORMAT(T 
GO  TO  20 
WRITE(6, 
FORMATU 
CONTINUE 
DO  201  I 
DO  201  J 
IF( IoLEo 
IFUoLEo 
IFUoGTo 
IFUoGTo 
WRITE(6, 
FORMAT (T 
GO  TO  20 
WRITEC6, 
FORMAT(T 
GO  TO  20 
WRITE(6, 
FORMAT(T 
GO  TO  20 
WRITEC6, 
FORMAT (T 
CONTINUE 
DO  202  I 
IFdoGTo 
WRITE(6T 
FORMAT(T 
GO  TO  20 
WRITE(6, 
FORMAT(T 
CONTINUE 


700) 

NAME'  ,T15t  'TEST*  ) 

701  ) 

ROWS' ) 

702) 

2,  «N      OBJF' ) 

=1,N 

9)     GO    TO    350 

703)    J 

2,'G       ACID 

0 

753)    J 

2,'G       A0,«,I2) 

=  1,M 

=  1»N 
90 ANDo JoLEo°)  GO  TO  351 
9oANDo  JoGT09)  GO  TO  352 
9oANDo  JoLEo9)  GO  TO  353 
9o ANDo JoGTo9)  GO  TO  354 
741)  ItJ 
2,  'L  A«-ti -•-, 
1 

742) 
2,  «L 
1 

743) 
2,  «L 
1 
744) 


tllf'.MD 


I»J 
A«  ,11,  S « ,12) 


I»J 
A'  ,12 


i»f »f  III 


ItJ 

A« 


tI2,«t  ",121 


=    ltM 
9)    GO    TO    320 
705)     I 

2t  'E      ARS  HI 
2 

72  0)     I 
2,»E       AR',I2) 


PUNCH    OUT    THE    COLUMNS    DATA    CARDS 
WRITE(6,706) 
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706      FORMAT! 'COLUMNS' ) 
JX    =    0 

JEND    =    M*N 

DO    300    J    =    1,JEND 

JX    =    JX    +    1 

00    302    K    =    1,M 

IF( JoLF0K*N,ANDo JoGTo (K-1)*N)     IX    =    K-l 

IFUoEQo  IX*N+1  )     JX    =    1 

302  CONTINUE 

DO    301    I    =    It ILIM 

IF(KEY(I,J)oNEol )    GO   TO    301 

M2    =    M    +    2 

DO    303    L    =    1,M2 

IF(  IoLEoL*NoANDo  IoGTo  (L-1)*N)     IA    =    L-l 

303  CONTINUE 

JA    =    I    -    IA*N 

IF(IoGTo(M+l)*N»     JA    =    I    -    (M+1)*N 

IF(  IXoLEo90ANDo  JX0LEo9oANDoI  A0  LE3Q,ANDo  JAo  LEo9)  GO  TO 
1360 

IF(IXoLEo9o  ANDo  JXoGTo9o  ANDo  IAo  LEo9,ANDo  JAo  LEo9)  GO  TO 
1361 

IF(IXoLEo9oANDoJX0GTo9oANDoIAoLEo9oANDoJAoGTo9)  GO  TO 
1362 

IF(  IXoLEo90ANDoJXoGTo9oANDo  IAoGTo9oAND3  JAo  LEo9)  GO  TO 
1363 

IF(  IXoLEo9oANDoJXoGT09oANDoIAoGTo9oANDo  JAoGTo9)  GO  TO 
1364 

IF(  IXoGTo90ANDoJXoGTo9o  ANDoIAoLEo9o  ANDo  JAoGTo9)  GO  TO 
1365 

IF( IXoGTo9oAND0 JXoGToQoANDoIAoGToQoANDo JAo LEo9)  GO  TO 
1366 

IF(IXoGT09oANDoJX0GTo9oANDoIAoGTo9oANDo  JAoGTo9)  GO  TO 
1367 

IF(  IXoLEo9oANDoJX0  LEo9oANDoI  Ao  LEo9oANDo  JAoGT09)  GO  TO 
1368 

IF( IXoGT09oANDo JXo LEo9o ANDoIAoGTo93 ANDo JAo LE,9)  GO  TO 
1369 

IF( IXoLEo9oANDoJX0LEo9oANDoIAoGTo9oANDo JAoGTo9)    GO    TO 

1370 

IF(lXoGTo9oANDoJX0LEo9oANDoIAoGTo9o  ANDo  JAoGTo9)  GO  TO 
1371 

IF(IXoLEo9oANDoJXoLEo9oANDoIAoGTo9o  ANDo  JAoLEo9)  GO  TO 
1372 

IF(IXoGT09oANDo  JX0LE09oANDoIAo  LEo9oANDo  JAoGTo9)  GO  TO 
1373 

IF( IXoGTo9oANDoJXoGT09oANDoIAoLEo9oANDo JAoLEo9)  GO  TO 
1374 

IF(  IXoGTo9oANDo  JXoLE09oANDoIAoLEo9oANDo  JAo  LEo9)  GO  TO 
1375 

360  IF(IoGTo(M+l)*N)    GO    TO    380 
WRITE(6,760)     I  X, JX , I  A, JA , AX(  I ,  J ) 

760  F0RMAT(T5,»X«  ,Il,,,,,Il,T15,,A«,Il»S'tIl,T2  5,F10o7) 
GO   TO    301 

380  WRITE(6,780)     IX, JX , JA , AX ( I » J ) 

780  FORMAT (T5,  ,X,,I1,,,,,U,T15,,AR',I1,T2  5,F1037) 
GO    TO    301 

361  IFUoGTo  (M  +  l  )*N)    GO    TO    381 
WRITE(6t761)     IX, JX, IA, JA,AX( I, J) 

761  FORMAT  (T5,  •X',Il,,,,,I2,T15,«A',Il,l,,,Il,T2  5,F10o7) 
GO    TO    301 

381  WRITE(6,781)  IX, JX , JA , AX ( I  ,  J  ) 

781  FORMAT (T5, •X',Il,*,,tI2,T15,,AR',Il,T25,F10o7) 
GO  TO  301 

362  IF(IoGTo(M+l)*Nl  GO  TO  382 
WRITE(6,762)  IX, JX ,  I  A, JA , AX( I , J ) 

76  2   FORMAT (T5, »X«,Il,l,«,I2,T15,,A',Il,«,»,I2,T2  5,F10o7) 
GO  TO  301 

382  WRITE(6,782)  I  X, JX , JA, AX ( I , J ) 

78  2   FORMAT (T5, «X«  ,I1,«  ,  • , I  2, Tl 5 , ' AR'  ,  I  2 , T25 , Fl Oo 7 ) 
GO  TO  301 

363  IF(IoGTo(M+l)*N)  GO  TO  383 
WRITE(6,763)  I  X, JX , I  A , JA , AX ( I , J ) 
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763 


383 
783 

364 

764 

384 
784 

365 

765 

385 
785 

366 

766 

386 
786 

367 

767 

387 
787 

368 

768 

388 
788 

369 

769 

389 
789 

370 

770 

390 
790 

371 

771 

391 
791 

372 

772 

392 
792 

373 

773 

F0RMATU5,  «X 
GO    TO    301 
WRITE(6,783) 
F0RMATU5,  'X 
GO    TO    301 
IFUoGT0(M  +  ! 
WRITE(6,764) 
F0RMAT(T5,« X 
GO    TO    301 
WRITE(6,784) 
F0RMAT(T5, ■ X 
GO    TO    301 
IFUoGTo  (M  +  l 
WRITE(6,765) 
F0RMAT(T5, «X 
GO    TO    301 
WRITE(6,785) 
F0RMAT(T5, «X 
GO    TO    301 
IF(  IoGTo(M  +  l 
WRITE(6,766) 
F0RMAT(T5, • X 
GO    TO    301 
WRITE(6,786) 
F0RMATU5, ■ X 
GO    TO    301 
IFUoGTo  (M  +  l 
WRITE(6,767) 
F0RMATU5,  'X 
GO    TO    301 

WRITE(6,787 
F0RMATU5,  «X 
GO    TO    301 
IFUoGT0(M+l 
WRITE(6,768) 
F0RMAT(T5,»X 
GO    TO    301 
WRITE(6,788) 
F0RMATU5,  '  X 
GO    TO    301 
IFUoGTo(M+l 
WRITE(6,769) 
F0RMATU5,  «X 
GO    TO    301 
WRITE(6,789) 
F0RMAT(T5,  «X 
GO    TO    301 
IF(IoGTo(M+l 
WRITE(6,77C) 
F0RMATU5, «X 
GO    TO    301 
WRITE(6,790) 
F0RMATU5, 'X 
GO    TO    301 
IFUoGTo  (M  +  l 
WRITE(6,771 ) 
F0RMATU5,  «  X 
GO    TO    301 
WRITE(6,791) 
F0RMATU5,  «X 
GO    TO    301 
IFUoGTo  (M  +  l 
WRITE(6,772) 
F0RMAT(T5, ' X 
GO    TO    301 
WRITE(6,792) 
F0RMATU5,  •  X 
GO    TO    301 
IF(  IoGT0(M+l 
WRITE(6,773) 
F0RMATU5,  fX 
GO    TO    301 


•til 

tf» 

IX, 

•t  11 

JX, 

IX, 
1  »I1 

GO 
JX, 

IX, 
'  »I1 

JX, 

)*N) 

IX, 

1  ,12 

GO 
JX, 

IX, 
1  »I2 

JX, 

)*N) 

IX, 

'  ,12 

GO 
JX, 

IX, 
'  f  12 

JX, 

)*N) 

IX, 

*  ,12 

GO 
JX, 

)  IX 
',  12 

,JX 

IX, 
1  ,11 

GO 
JX, 

t1  f 

IX, 
•  1 11 

JX, 

)*N) 

IX, 

«,I2 

GO 
JX, 

IX, 
•♦12 

JX, 

)*N) 

•Jfi 

GO 
JX, 

IX, 
*  t  11 

JX, 

)*N) 

IX, 

1  tI2 

GO 
JX, 
»•  t 

IX, 
1  »I2 

JX, 

»*N) 

IX, 

•til 

GO 
JX, 

>\h 

JX, 

)*N) 

IX, 

'  f  12 

GO 
JX, 

JA, AX(I  ,J 
SI2,T15, 

TO    384 
IA, JA, AX( 
f,I2,T15, 

JA,  AXU  ,J 
SI2,T15, 

TO    385 
IA, JA,AX( 
M2.T15, 

JA,  AXU  ,J 
SI2,T15, 

TO    386 
IA, JA,AX( 
•tI2,T15, 

JA,  AXU  ,J 
SI2,T15, 

TO    387 
IA, JA,AX( 
S I2tT15t 

, JA,AX( I  , 
•»I2,T15, 

TO    388 
IA, JA, AX( 
SI1,T15, 

JA,  AXU  ,J 
•tlltTiS, 

TO  389 
IA, JA,AX( 
•tIltT15i 

JA,  AXU  ,J 
SI1,T15, 

TO    390 
IA, JA,AX( 
•fIliT15i 

JA,  AXU  ,J 
MltT15, 

TO    391 
IA, JA,AX( 
'»I1»T15, 

JA,  AXU  ,J 
•tIl»T15, 

TO    392 
IA, JA,AX( 
•»I1»T15, 

JA,  AXU  ,  J 
•f IlfTlSi 

TO    393 
IA, JA, AX( 
SIltT15, 


AR 


,  J 

A« 


AR 


t  J 
A« 


AR 


,  J 

A' 


AR 


,  J 

A' 

) 

AR 


,  J 

A« 


AR 


tJ 
A' 


AR 


.  J 
A* 


AR 


,  J 

A" 


AP 


t  J 

A' 


AR 


,  J 

A» 


12, 


til 


12 


,12 


lit 


i  12 


12, 


,11 


12, 


,12 


Il» 


iI2 


12 


til 


12, 


,12 


12 


,12 


12, 


,11 


,• ,Il,T25,F10o7) 


T25,F10o7) 


,• ,I2,T25,F10o7) 


T25,F10,7) 


, »,I2,T25,F10o7) 


T25,F10o7) 


, ' , Il,T25,F10o7) 


T25,F10o7) 


, ' ,I2,T25,F10o7) 


T25,F1Cd7) 


, «, I2,T25,F10o7) 


T25,F10o7) 


,  SIl,T25,F10o7) 


T25,F10o7) 


,,,I2,T25tF10o7) 


T25,F10o7) 


, ',I2,T25,F10o7) 


T25,F10o7) 


,« ,Il,T25,F10o7) 


T25,F10o7) 


II,' , «,I2,T25iF10o7) 


44 


393  WRITE(6,793)     I  X, JX , JA, AX ( I , J ) 

79  3      FORMAT  (TS^XS  12,'  ,  ' ,  1 1  tU5,  •  AR«  ,  I  2  tT25,  FlOo  7) 
GO    TO    301 

374  IF( IoGT0(M+l )*N)    GO    TO    394 
WRITE(6,774)     I  X,  JX  ,  I  A,  JA  ,  AX  (  I  ,  J  ) 

774      FORMAT (T5t • X« , 12 , ■ , • , I 2,T1 5 , ' A« , II , • , • , II , T25,F10o7) 
GO    TO    301 

394  WRITE(6,794)  I  X,  JX  ,  JA  ,  AX  (  I  ,  J  ) 

794   FORMAT (T5, ,X',I2,',,,I2,T15,,AR«,Il,T?5,F10o7) 
GO  TO  301 

375  IF(IoGTo(M+l)*N)  GO  TO  395 
WRITE(6,775)  IX, JX , I  A, JA , AX( I , J ) 

77  5   F0RMAT(T5,'X» , 12 , ■ , « , 1 1,T15 , ■ A' , II , ■ , « , II ,T25,  FlOo 71 
GO  TO  301 

395  WRITE(6,795)  IX, JX , JA , AX( I , J ) 

79  5   FORMAT (T5, ,X«,I2,,,«,Il,T15,'AR',Il»T25,F10o7) 
301   CONTINUE 

IF( JoLToJEND)  GO  TO  300 

DO  305  IZ  =  1,M 

JZ  ■  IZ 

IF(IZoGTo9)    GO    TO    304 

WRITE(6,709)     IZ,     JZ 
709       FORMAT (T5, •  XA«  ,I1,T15,  'OBJF* ,T25,< la  •  ,T40, »AR« ,I1,T50, 

l'lo'   ) 

GO    TO    305 

304  WRITE(6,708)     IZ,    JZ 

70  8       FORMAT (T5,  •XA,,I2,T15,,0BJF»,T25,,lo' ,T40,  «AR»,I2,T50, 
l'lo'  J 

305  CONTINUE 
300      CONTINUE 

PUNCH  OUT  THE  RIGHT  HAND  SIDE  DATA  CARDS 


R 
R 
R 
R 
R 
P 
R 
R 


710 


717 

333 
718 
325 


311 
711 

312 
722 

313 
72  3 

314 
724 
30  6 


712 

315 
72  5 
30  7 

713 

801 


WRITE(6, 
FORMAT( • 
DO  325  I 
IFdoGTo 
WRITE(6, 
FORMAT(T 
GO  TO  32 
WRITE(6, 
FORMAT<T 
CONTINUE 
DO  306  I 
IF(MoLEo 
IF(M0LEo 
IF(MoGTo 
IF(MoGT0 
WRITE(6, 
FORMAT(T 
GO  TO  30 
WRITE(6, 
FORMAT (T 
GO  TO  30 
WRITE(6, 
FORMAT(T 
GO  TO  30 
WRITE(6, 
FORMAT(T 
CONTINUE 
DO  307  I 
IFdoGTo 
WRITE(6, 
FORMAT(T 
GO  TO  30 
WRITE(6, 
FORMAT(T 
CONTINUE 
WRITE(6, 
FORMAT ( ■ 
IF( ITEST 


710) 
RHS'I 

=    1,N 
9)    GO    TO   333 

717)  I,     PX(I,1) 

5, •Bl« ,T15, 'AO, ' ,Il,T25,F10o7) 
5 

718)  I,    PX(I,1) 

5, •El« ,T15, «A0,« ,I2,T25,F10o7) 

=    1,N 

9oANDoIoLEo9)     GO    TO  311 

9oANDoIoGT09)     GO    TO  312 

9oANDoIoLEo9)     GO    TO  313 

9oANDoIoGT<>9)    GO    TO  314 

711 )  M,     I,     X(I  ,1) 

5,  'El' ,T15, «A« ,11, ',', Il,T25,F10o7) 
6 

722)  M,     I,     X(I,1) 

5,,Bl«,T15,«A»,Il,«,«,I2,T25,F10o7) 
6 

723)  M,     I,     X(I,1) 

5,«Bl,,T15,»A',I2,,,',Il,T25,F10o7) 
6 

724 )  M,     I,     X(I  ,1) 
5,'Bl,,T15,«A',I2,,,,,I2,T25,F10o7) 

=    1,M 
Q)    GO    TO   315 

712)  I 

5,  «B1*  ,T15,  •ARSI1.T25,  '  lo  •  ) 
7 

725)  I 

5, «E1« ,T15, «AR» , I2,T25,  •  lo  f  ) 

71?) 

ENC£TA« ) 

oEQoIRl )  GO  TO  999 
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PUNCH    OUT    THE    APPROPRIATE    RHS    REVISION    DATA    CARDS 


714 

308 
715 
309 

716 


R 
R 
R 
R 
R 
R 


334 
326 


321 

322 

323 

324 
310 

999 


IF( ITEST 

WRITE(6, 

FORMAT( • 

GO    TO    30 

WRITE(6, 

FORMAT( » 

WRITE(6, 

WRITE(6, 

FORMAT(T 

IT1    =    IT 

DO    326    I 

IFUoGTo 

WRITE(6t 

GO    TO    32 

WRITE(6, 

CONTINUE 

DO    310    I 

IF(MoLEo 

IF(M0LEo 

IF<MoGTo 

IF(M0GTo 

WRITE(6, 

GO    TO 

WRITE( 

GO    TO 

WRITE( 

GO    TO 

WRITE( 

CONTINUE 

WRITE(6t 

ITEST    = 

GO    TO    80 

CONTINUE 

STOP 

END 


31 
6, 
31 
6, 
31 
6, 


oGTo9)    GO    TO    308 

714)  ITEST 

NAME1  ,T15t  'REV  ,11  ) 
Q 

715)  ITEST 

NAME'  »T15,  'REV  »  12) 

71 C) 

716) 

3,  "MODIFY*  ) 

EST    +    1 

=    1,N 
9)    GO    TO    334 

717)  I,PX(I,IT1) 
6 

718)  I,     PXU,IT1) 

=    1,N 

9oANDo IoLEo9)    GO    TO  321 

9oANDo  IoGT09)     GO    TO  322 

9o ANDo IoLEo9)     GO    TO  323 

9o  ANDo  IoGTo9)     GO    i"o  324 
711)     M,     I»     X(ITIT1) 


0 

722) 

0 

723) 

0 

724)    M, 


M, 


It 
It 

It 


X(I ,IT1  ) 
X(I ,IT1) 
X(I ,IT1) 


713) 
ITEST    + 
1 


^^ij:^!!!)):^^*^^**'''********5!1******* 
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****************************** 
*  *  DATA  GENERATOR  -  TEST  FOR  A  REACHABLE  SET  IN  M  STEPS  *  * 
******  ROBERT  Lo  ARMACOST  -  -  -  1  JULY  1970  ****** 
*************  ***************** 

****************************** 

THIS  PROGRAM  WILL  PUNCH  OUT  A  SET  OF  DATA  CARDS  FOR  A 
LINEAR  PROGRAM  USING  THE  MPS  PACKAGE  ON  THE  IBM  360/67, 
THE  BASIC  INPUT  TO  THE  PROGRAM  IS  A  SQUARE  PROMOTION 
MATRIX  OF  SIZE  N,  AND  K  VECTORS  WHICH  ARE  TO  BE  TESTED 
TO  DETERMINE  IF  THEY  ARE   ELEMENTS  OF  THE  SET  B<M>o 
THE  PROGRAM  WILL  GENERATE  DATA  DECKS  FOR  A  SERIES  OF 
REACHABLE  SETSo   MBEG  IS  THE  STARTING  INDEX  FOR  THE 
SUBSCRIPT  OF  THE  TESTo   MEND  IS  THE  ENDING  INDEX  FOR 
THE  SUBSCRIPT  OF  THE  TESTo   THESE  INDICES  ARE  READ  IN 
AS  DATAo   THE  FIRST  VECTOR  TO  BE  TESTED  IS  CALLED 
"TEST"o   EACH  SUCCEEDING  VECTOR  IS  CALLED  "REVI"  WHERE 
I  IS  THE  NUMBER  OF  THE  REVISION  OF  THE  RIGHT  HAND  SIDEo 
THE  VARIABLE  IR  IS  THE  NUMBER  OF  REVISIONS  IN  THE 
CONTROL  PROGRAMo   THE  ARRAYS  SHOULD  BE  DIMENSIONED  IN 
THE  FOLLOWING  WAY  - 

AX((MEND)*(N+11  +  N,  MENDMN  +  1)) 

KEY  -  S£ME  AS  AX 

X(N, IR+1) 

P(N,N) 

PT(NtN) 

PTI(N»N)o 

THE  PROGRAM  WILL  ACCOMODATE  A  PROMOTION  MATRIX  WITH 
N  =  99  AND  WITH  MEND  <  100o 

****************************** 

****************************** 

THIS  PROGRAM  IS  CONSTRUCTED  BY  REMOVING  THOSE  CARDS 

MARKED  WITH  AN  R  IN  PROGRAM  TO  TEST  FOR  RETURN  IN  M 

STEPS,  AND  TC  REPLACE  THOSE  MODULES  MARKED  WITH  NUMBERS 

BY  THE  MODULES  LISTED  BELOW  WITH  THE  CORRESPONDING 
NUMBERSo 

****************************** 

1  READ(5,600)  MBEGt  MEND,  N,  IR 

2  CALCULATE  THE  TRANSPOSE  OF  THE  PROMOTION  MATRIX,  PT, 
2     AND  CALCULATE  I  -  PT  WHERE  I  IS  THE  IDENTITY  MATRIX 
2 

2  DO  20  1=1, N 

2  DO  20  J=1,N 

2  •  PT(I  ,J)  =  P(J,I  ) 

2  PTKI.J)  =  -PT(  I  ,J  ) 

2  IF(IoEOoJ)  PTI(I,J)  =  lo  +  PTI(I,J) 

2  20  CONTINUE 

3  INSERT  PTI  IN  THE  FIRST  N  ROWS  AND  COLUMNS  OF  THE 
3     CONSTRAINT  MATRIX,  AX 

3 
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3  DO    30    1=1, N 

3  DO    30    J=ltN 

3  AXU  ,J)     =    PTK  I,  J) 

3  30    IF(AX( I, J)oNEoO)    KEY(I,J) 

4  (M-1)*N    ROWS    OF    AX 


=    1 


$#$'-jc$:£$$$)t';:4c^c$$$$£$$£'>fc3t':#)j($$$$$4: 
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****************************** 

*  STEADY  STATE  SET  -  EXTREME  POINTS  AND  TEST  FOR  INCLUSION  * 

******  ROBERT  Lo  ARMACOST  -  -  -  1  JULY  1970  ****** 
****************************** 

THE  MAIN  PROGRAM  IS  DIMENSIONED  TO  HANDLE  A  PROMOTION 
MATRIX  WITH  15  STATES,   THE  PROMCTI ON  MATR I  X  IS  P  AND 
ITS  INVERSE  WHICH  HAS  BEEN  NORMALIZED  IS  PINo   THE 
ROWS  OF  PIN  ARE  THE  EXTREME  POINTS  OF  THE  STEADY  STATE 
SETo 

THE  EXTREME  POINTS  OF  THE  STEADY  STATE  SET  ARE  USED 
TO  DETERMINE  IF  A  TEST  POINT  IS  IN  THE  STEADY  STATE 
SET  BY  CALCULATING  THE  NECESSARY  COEFFICIENTS  TQ 
MAKE  THE  TEST  POINT  A  LINEAR  COMBINATION  OF  THE 
EXTREME  POINTSo   IF  THE  SUM  OF  THE  COEFFICIENTS  IS 
GREATER  THAN  ONE,  OR  IF  ANY  COEFFICIENT  IS  NEGATIVE, 
WE  CONCLUDED  THAT  THE  TEST  POINT  IS  NOT  IN  THE  STEADY 
STATE  SETo   THE  COEFFICIENTS  ARE  DENOTED  ALF,  AND  THE 
TEST  POINTS  APE  X,   THE  PROGRAM  IS  SET  UP  TO  ACCEPT, 
FIVE  SUCH  TEST  VECTORS  WITH  FIFTEEN  RANKSo 

THE  SIZE  OF  THE  PROMOTION  MATRIX  AND  THE  LENGTH  OF  THE 

TEST  VECTOR  ARE  INPUT  TO  THE  PROGRAM  IN  THE  VARIABLE 

NNo   THE  NUMBER  OF  TEST  VECTORS  IS  INPUT  IN  THE 
VARIABLE  KKo 

*****  ************************* 


DIMENSION  UNI (15,15) ,P(15, 15),  P 
1M(15),  SUM(15),  ALFU5),  X(15,5) 
READ(5,599)  NN,  KK 
599   FORMAT(2I5) 

CALL  WORK(UNI,P,PIN,L,M,ALF,X,NN,KK,SUM) 
CALL  EXIT 


IN(15,15),     LQ5) 


END 


****************************** 

*  *    *    SUBROUTINE    WORK    -    THIS    SUBROUTINE    CALCULATES    THE    *    *    * 

*  *    *    EXTREME    POINTS    OF    THE    STEADY    STATE    SET,     M,     AND         *    *    * 

*  *    *    TESTS    EACH    TEST    VECTOR    FOR    INCLUSION    IN    THAT    SET    *    *    * 


SUBROUTINE    WORM  UN  I , P , PIN , L , M, ALF, X ,NN , KK, SUM) 

DIMENSION    UM (MN,NN) ,P(NN,NN) , PIN(NN,NN) ,L(NN), 
1M(NN),     SUM(NN),     AL F ( NN ) , X( NN , KK) 

RE AD (5, 600)     <(P(I,J),J=1,NN),I=1,NN) 
6C0       FORMAT(8F10oO) 

RE AD (5, 601 )     ( ( X( J,K), J=1,NN) ,K=1,KK) 
601       FORMAT(8F10oO) 

INITIALIZE 

DO    10    I    =    1,NN 
L(I)    =    0 
M(I)    =    0 
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SUM( I)  =  Oo 
DO  10  J  =  1,NN 
UNK  It  J)  =  -  P(I  TJ) 

IF(IoE0oJ)  UNIU.J)  =  lo  ♦  UNI(ItJ) 
10    CONTINUE 

WRITE(6,699) 

WRITE  OUT  PROMOTION  MATRIX 

DO  2  1=1, NN 

WRITE(6,700)  (P< I, J) ,J=1,NN) 
WRITE(6, 710) 
2     CONTINUE 

CALCULATE  THE  INVERSE  OF  P 

CALL  MINV(UM  ,NN,D,L,M) 

CALCULATE  THE  NORMALIZED  INVERSE  OF  P,  PIN 

DO  20  1=1, NN 

DO  15  J=1,NN 

SUM(  I)  =  SUM(  I)  +  UNK  I,  J) 
15    CONTINUE 

DO  25  J=1,NN 

PIN(  I,  J)    =    UNKI  ,J)    /    SUM(  I) 
25         CONTINUE 
20         CONTINUE 

WRITE(6,699) 

699  F0RMATUH1) 

WRITE  OUT  NORMALIZED  INVERSE  OF  P 

DO  5  1=1, NN 

WRITE(6,700)  (PINU,  J)  ,J  =  1  ,NN) 
WRITE(6,710) 
5     CONTINUE 

700  F0RMAT(15F8o6) 
710      FORMATUH    ) 

TEST  TO  SEE  IF  TEST  VECTOR  IS  IN  THE  STEADY  STATE  SET 

DO  130  K=1,KK 

MM  =  K 

DO  1  1=1, NN 
1     ALF( I)  =  Oo 

ALF(l)  =  X(1,K)  /  PIN(1,1) 

TEST  =  ALF(l) 

JEND  =  NN-1 

DO  120  J=2,JEND 

IEND  =  J  -  1 

SUMD  =  Oo 

DO  110  1=1, IEND 
110   SUMD  =  SUMD  +  ALF(I)  *  PINU, J) 

ALF(J)  =  (X(J,K)  -  SUMD)  /  PINU, J) 

TEST  =  TEST  +  ALF(J) 

IF(ALF(J)oLTo0)  GO  TO  100 

IF(TESToGTol)  GO  TO  200 
120   CONTINUE 

ALF(NN)  =  lo  -  TEST 

IF(ALF(NN)oLTo0)  GO  TO  100 

IF(ALF(NN)oLT0X(NN,K) )  GO  TO  200 
130   CONTINUE 

WRITE(6,702) 

702  F0RMAT(1H1,T20,'THE  VECTOR  IS  IN  THE  STEADY  STATE  SET* 
1//T2C,'ALF  COEFFICIENTS' ,T45, 'TEST  VECTOR'//) 

DO  201  LL=1,NN 

WRITE(6,703)  ALF(LL),  X(LL,MM) 

703  F0RMAT(T20,F10o8,T45,F10o8) 
201      CONTINUE 

IF(KoLToKK)    GO    TO    130 
GO    TO    999 


50 


100       WRITE<6,704) 

704  FORMATUH1  ,T20t  'THE    VECTOR     IS 
It*     SET,     ALF    IS    NEGATIVE'//T20, 
2'TEST    VECTOR'//) 

DO    202    LL=1,NN 

WRITE<6,703)     ALF(LL),     X(LL,MM) 

202  CONTINUE 
IF(KoLToKK)    GO    TO    130 
GO    TO    999 

200       WRITE(6,705) 

705  FORMAT  JdHl,T20»«  THE    VECTOR    IS    NOT    IN    THE 
1,'     SET,     TEST    IS    GREATER    THAN    0NEV/T20, 
2'ALF    COEFFICIENTS' ,T45, 'TEST    VECTOR'//) 

DO    203    LL=1,NN 

WRITE(6,703)    ALF(LL),     X(LL,MM) 

203  CONTINUE 
IF(KoLToKK)    GO    TO    130 

999   RETURN 
END 


NOT  IN  THE  STEADY  STATE' 
•ALF  COEFFICIENTS' »T45, 


STEADY  STATE' 


*^c**s{c3{t********jJc**^:A***3jc******* 
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